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Secular Perturbations of the Planets. 

Br G. W. Hill. 



Gauss first clearly indicated the role elliptic functions play in this subject.* 
Halphen has since presented the investigation in a very elegant manner, f The 
modifications made by the latter in the procedure of Gauss are chiefly the trans- 
ference of the origin of rectangular coordinates to the attracted planet, and, instead 
of the differential of the eccentric anomaly, the adoption of the element of area 
described by the radius of the disturbing planet expressed in terms of the differ- 
entials of the rectangular coordinates. He also appeals to the qualities of the 
cone formed by the orbit of the attracting planet as contour of base and the 
position of the attracted planet as vertex ; this improvement, however, had been 
previously indicated by Bour. J A remarkable degree of elegance is attained by 
these changes; but it seems to me that additional statements are needed to 
show the connection with the astronomical problem which originally suggested 
the investigation; for Halphen, like Gauss, treats only the attraction of a certain 
form of ring. This, of course, is to ignore the second integration which the 
problem demands. Perhaps, therefore, I shall be pardoned if I here attempt to 
supply the mentioned lack. 

In the fashion of Halphen we take the attracted planet as the origin of rectan- 
gular coordinates, but the orientation of the axes is, for the present, left indeter- 
minate. The coordinates of the attracting planet we denote by x , y , z ; and the 
coordinates of the Sun, which are the negatives of those of the attracted planet 
referred to the Sun, will be x , y , z . Let p be the distance of the attracting 
planet from the origin, so that f = a? + y* + 2 s ; and let g denote the planet's 

* Gauss, Werke, vol. Ill, pp. 331-355. 

f Q. H. Halphen, Trait6 des Fonctions Elliptiques efc de leurs Applications, Tom. II, pp. 310-328. 

X Journal de l'^cole Polytechnique, Cahier XXXVI, pp. 59-84. 
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mean anomaly. Then the secular perturbations of the attracted planet depend 
on the three definite integrals 

j-r*ag, j-pJLm, ±r^d g . 

But while g is passing from to In, the area described by the radius of the planet 
augments from to nab, if a and b are severally the major and minor semi-axes. 
Thus, if a denote this varying area, the preceding integrals may be written 



ndb y 



nab Jo p 3 ' nab Jo p 3 ' ' nab J° p 3 



da. 



The tetrahedron, with da as base and the origin as vertex, has, for volum e 
one sixth of the following expression : 

6 V = x (ydz — zdy) + y (zdx — xdz) + z (xdy — ydx). 

But, if h denote the perpendicular from the origin on the plane of the orbit of 
the attracting planet, we also have 3 V = hdo. Hence 

do = ll ( ydz ~~ zdy) + ^ { - zdx ~ xdz) + "S" ( xdy ~ ydx) • 

Then, if we derive the quantities P x , P y ,P t , etc., from integrating the expressions 
j d -i ^jydz — zdy) , p _ x y{ydz- zdy) , p _ , zjydz-~_zdy) 

9 r r 

, rt „.x (zda — adz) , n _,J (z^a — crdz) , n _, z (zcfa; — a-cfe) 

« V* — $ ts ' a "y ~ 2 rs ' a V* — * ^3 ' 

dK=i ^dy^ydx) j dRv=l y{xdy-ydx) ) ^ = ^ zjxdy-ydx) ^ 

9 9 y 

around the whole orbit, the integrals above, which we will denote by X, Y, Z 
will be given by the following expressions: 

X = ~nabh ^° P * + y *®' + Z ° R ^ ' 

Z = ^ (z P, + y & + z #*) . 

The nine quantities involved in these expressions and obtained through integration 
are homogeneous and of the dimension zero with respect to the linear unit. 
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Moreover, if, in the cone formed by the orbit of the attracting planet as directrix 
and the origin as vertex, the plane of the base is shifted in any manner whatever, 
these nine quantities remain unchanged. For example, consider the expressions 



x (ydz — zdy) y (ydz — zdy) z (ydz — zdy) 
If we put 



P 3 "~" ' ~ P 3 ' ' "" P 3 



x = p sin 6 , y = p cos cos A. , z = p cos 6 sin X , 

they are transformed into 

sin 6 cos 2 9 dk , cos 3 6 cos A d A , cos 3 sin a d a 

Let the equation of the cone be 

Ax* + By* + Cs 2 f Ztyz + Fzx + Fxy = 0, 
or 

J.tan 2 + J? cos 2 a + C sin 2 a, + D sin a, cos a, + Et&nd sin a + i^tanflcos a =0. 
Since, from this equation, is obtainable as a function of a, it is evident that the 
integrals of the preceding differential expressions, extended to the whole course 
of variation of X, depend solely on the elements of the cone and are altogether 
independent of the plane section called the base. 

A simple addition of the differentials shows that 

P x + Q y + R t = 0. 
Also we have 

d(Q x -R y ) = \d(^-), d{R x -P,) = \dQL), d(P y -Q x ) = hd(^-y 

But the second members of these equations, integrated along the orbit to the 
point of beginning, give zero as the result : hence 

Thus 

X = nabh ^ X ° Px + y °® x + Z ° P ^ ' 

Y ~ TtoM I**®' + y °® y + Z ° Ry ) ' 

Z ~ nabh ^ Pz + y ° By + Z ° R ^' 
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And, if we put 

* = 2^6A tx ° p * + ylQ y + Z *° R ' + 2x <v° Q ' + 2mR y + 2Zo!BoP ^ ' 

we shall have 

X= §*- F= — Z = — 

8a5 ' 3y ' dz ' 

The orientation of the axes of coordinates which serve to define the variables 
x, y, z, x , y , z has been left undetermined ; but now suppose that the axes of 
symmetry of the cone are employed for this purpose. Then the equation of the 
cone takes the form 

*L + jL + A =o 

<X % (X y (X e 

G x , G y , G z being constants of which two are of one sign and the other of the 
opposite sign. Plainly, if this equation is satisfied by the set of values x, y, z, 
it is satisfied by any of the eight sets ± x, ±y, ±z. Consequently, each posi- 
tive element of the six quantities Q z , B y , B x , P z , P y , Q x is accompanied by a 
corresponding negative element. Thus, in this case, these quantities vanish. 
"With this selection of axes we, therefore, have 

X=-%-P x , Y=-^rQy, Z=-\ l B z . 
■jtabh °° 7tabh^ y nabh * 

The naming of the coordinates is, of course, arbitrary, but, to settle the 
choice, we suppose that G x , G y , G z are in the order of algebraic magnitude, the 
first and second being negative, while the last is positive. The equation of the 
cone appears to involve three variables, but, as we may divide the left member 
by the square of any one of them, it is, in reality, a relation between two vari- 
ables ; thus, but one variable can be regarded as independent. The equation is 
then satisfied if we make 



a; = 6*/— G x cos T, y = eV— G v sin T, z = e«/G z> 

where e is an indeterminate which disappears when the substitution is made in 
dP x> dQ v , dB 2 , and T is the new variable introduced by Gauss and may be 
regarded as indicating the position of the planet in its orbit ; its function, in this 
respect, being precisely similar to those fulfilled by the mean, eccentric and true 
anomalies, and it may thus, with propriety, be designated as a perspective anomaly. 
When 7 goes from to In, the planet makes a complete circuit of its orbit. 
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The substitution made, we have 



ydz — zdy = — e 2 V — G y G z cos TdT, 
zdx — xdz = — e 2 */ — G x G z sin TdT, 
xdy — ydx = e 2 V G t G y d T, 

p 2 = e % [G e — G y sin 2 T— G x cos 2 T] . 

The quantities G x , G y , G z are usually determined in such a way that G K G y G z 
= atWih also we may introduce k the modulus of the elliptic integrals involved 
and m such that 



*ir 



Then our integrals take the forms 

x _ ^^_l /*■! cos 8 y^y y _ yo_ _2 /*.! sin 2 y#y 

m 8 nJ° (l-Fsin 2 ? 7 )!' m 3 nJ° (l—7^mi 2 T)V 

<m i ft Jo (1 — , 



(1 — & 2 sin 2 y)|* 

The methods of evaluating these definite integrals severally proposed by 
Legendre, Gauss and Jacobi have all about the same degree of rapid convergence 
but that of the last is to be preferred because it expresses the values explicitly in 
terms of a parameter q called the nome. Putting h = sin $, q can be derived 

from the equation 

q + q 9 + <f 5 + • • • • r sin J -f 

l + 2(2* + 2 16 +2 36 + ....) Ii + V~cos0 J' 

The solution is most readily accomplished by the method of tentation. Then, if 
we adopt two functions of q, K and L such that 

f __ (1 + V cos ft) 3 g— 4g 4 +9g 9 — 16<? 1,! + ••■• 
sin 2 ftcosfft [l+2( ? * + ? 16 + ••••)] 3 ' 
we have 

tn 3 in 3 v ' m 3 v ' 

It will be seen from these expressions that, if X, Y, Z are regarded as the 

* For the proof of these formulas, reference may be made to Bertrand, Calcul Integral, Liv. III., 
Chap. VII. 
42 
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components along the axes of coordinates of a force acting on the attracted planet, 
one elliptic integral suffices for determining the orientation of the resultant, but 
that an additional one is required if the magnitude of the latter is to be found. It 
will be an advantage, therefore, if instead of tabulating if and L as functions of 
h, q or 6, we take two other quantities J!f and x, such that 

M= K — L sir 0, sin 2 x = -^ T . i6 , cos 2 x - 



/f— Zsin 2 0' K—Lsin*0 

Then we shall have 

JC = 5-surx. x , Y= — — o-cos'jc. ?/ , Z — — „ z . 

m 6 mr ^ m 3 

Let B denote the magnitude of the resultant, H and A severally the latitude 
and longitude of the point in the heavens towards which it is directed ; the circles 
of reference being the principal axes of the sphero-conic traced in the heavens 
by the frequently mentioned cone. Also let r denote the distance of the Sun 
from the attracted planet and yi , a, a severally its latitude and longitude referred 
to the same circles. Then our equations will stand 

M 

R cos H sin A = 3- sin 2 x . r cos yj sin X, , 

It sin IT — — — T cos 2 x . r sin v> , 

rnr 



R cos H cos A = — * . r cos >? cos a, a . 



m° 






If we put 



we shall have 

, . • 2 . ~ . tt cos 3 x tan Yi a cos A D N cos \ 
tan A = — sin^ x tan ^, , tan H = i5 , R = = — 2_ . 

cos a cos H cos A 

As, except in very particular cases, it is not easy to select at the outset the 
axes of symmetry of the cone for the adopted axes of coordinates, we must find 
the position of these axes in reference to another system which is known. The 
equation of the cone having a very complicated expression when the axes of 
coordinates are quite general, we select a particular system such that the treat- 
ment may be as easy as possible. Let the axis of x have the orientation of the 
line going from the centre of the ellipse described by the attracting planet to its 
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perihelion, that of?/ the orientation of the line going from the centre to the point 

where the eccentric anomaly is 90° and that of z the orientation of the line going 

from the centre to the north pole of the plane of the orbit. Let the coordinates 

of the centre of the ellipse be, in their order, A, B, G; we prefer these designations 

although, with e as the eccentricity and the previous notation, they have the 

equivalents 

A — ;r — ae, B = y , G = z . 

Then the equations of the orbit of the attracting planet are 

A\* , fy~B\ 2 _ 



(^)+(^)=>' ■=" 







a ~y 


— 


B ~] 


g ? 


X - 


-~G Z 




y- 


~-G z 


G* 


"*. 


a 




— 


b 



In order to have the equation of the cone so frequently mentioned it is only nec- 
essary to multiply the several terms of the first equation by factors selected from 

the equivalent quantities 1 = ~ = -^ , in such a way that they may all become 

homogeneous and of two dimensions in x, y, z. Thus the equation of the cone 
may be written in the shape 



0. 



This equation is not, iu general, referred to the axes of symmetry, hence we 
proceed to make the linear transformation of variables which will bring this 
about. Let x, T, z denote the rectangular coordinates referred to the sym- 
metric axes of the cone, and write the formulas of transformation thus : 

x = ax + Py + yz , a? -f 0» + f = 1 , aa' + pP' + yy 1 = 0, 
Y = a'x + P'y + y'z , a' 8 + P P +y" = 1, a'a" + P'p" + //' = 0, 
z = a"x + P"y + y"z , a" 9 + p"* + y»' = 1 , a"a + P"P + y"y = , 

the inverse of which are 

a: = aX + a'Y + a"z, a 2 + a' 2 + a" 2 = 1 , a/3 + offi + a"P" = , 
y = pX + P'Y + P"Z, P* + P<* + P»* =z 1 , Py + py + P"y" = , 
z=yX + y'Y + y"z, y 2 + y* + y"* = 1 , y a + y'a' + y"a" = 0. 

The equation of the cone, after substitution of the new variables, should 
take the form 

X 2 . Y 2 , Z 



G-% 



+ 



G„ 



+ ^=o, 
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where G x , G y , G z are constants and functions of the quantities A, B, G , a,b. 
They are the roots of a certain cubic which may be obtained in the following 
way : Let V= denote the first form of the equation of the cone, and set 

_^_~i<l? jl — i <*r z —i 37 

G _t dx ' G ~^ dy ' G ~* dz ' 

From these equations, which are linear in x, y, z, eliminate these variables; the 
result is a cubic in G whose roots are the values of G x , G y , G z . In the special 
form of V with which we have to deal, these equations, after a slight modifica- 
tion, can be thus stated : 

x jl z y B z z z __ A x By 

G "" ~C G + a* ' ~G~~G G+b* ' ~G ~~ W ~0 ~G ~ ~G ~G ' 

The elimination ofx,y,z from these equations gives 

A* B* ,&_ 

,2+ a _i_ hi + a. * • 



G + a?^G+b*^G 

Tt is known that the roots of this equation in G are all real, two of them being 
negative and one positive. The mode of assignment of these roots as values of 
G x , G y , G, has already been described. Thus we are enabled to discover the 
values of the latter without reference to the nine coefficients a, (3, y, etc. But in 
our further progress we shall need to know the latter. They are readily found 
from the first and second equations of the penultimate group combined with the 
conditions a 2 -f /3 3 + y a = 1, etc. In the first we make x = a, y = (3, z= y and 
set Q x for G ; again make x =■ a', y = /?', z = y' and set G y for G ; lastly make 
x =. a", y = /3", z = y" and set G z for G. Thus we have the equations 

A G x n B G x 

a ~-G G^+^ r ' P ~ ~C G x + b* y ' 



"*" C* KG x +a*J h C* \G~+bV 
~0 G y + a* / ' y ~ ~0 G y + b 



I -a. G v , ni B 6r„ 



y>* 



<x!>: 



i i- At ( G « Y i B V G * Y ' 

— ®J. n ," Rii — & ®* « I' 



+ aV ' C 2 \G y + V. 

--. - v " &'—— ""* y 



y" 2 : 



+ (7* U, + a* J ^ G* \G~+VJ 
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These equations furnish the values of the nine coefficients of the transformation. 
The signs of the three y's are indeterminate : one may take them as positive. 
The cubic in G furnishes us with the relations 

' n i 1.2 i n — L > 



G K + a* ' G x + b* ' G x 

~ n i 2.2 ~ si x > 



G y +a? ' G y + b* ' G y 

A 2 B 2 W _ 

i si i i2 ' n * ' 



#, + a 2 ^ 0, + 2> 2 ' ^ 

Eegarding J. 2 , -B 3 , C 2 as the unknowns in these equations, their solution gives 

4 s_ (£« + <)(ff j , + a 2 )(fl, + a 2 ) 
a 2 (a 3 — J 2 ) 

B B - (ff, + ft , )(ffl, + s |, )(g, + &») 

~~ 6 2 (6 2 — a 2 ) 

/72 __ Cfj, ur y Cr z 



2 7,2 



a 2 # 

By means of these values we can eliminate J., B, G from the foregoing expressions 
for the coefficients of transformation, and thus obtain values depending only on 
the five quantities G x , G v , G z , a, b. We have thus 



a 



G x {G x + V)(G v + a?){G z + a*) 



a?(a*-b*) {G x - G y ){G x -G z ) 



R , _ 1 G x (G x + a*)(G y + b*)(G z + b») 

P -6»(6»_o»)~ { G x -G y ){G x -G z ) 
-_ 1 G y G z {G x + a?){G x + V) 
r ~aW {G x ~G y ){G x -G z ) • 

The values of a' 2 , /3' 2 , y' 3 and again of a" 3 , /3" 3 , y" % , are obtained from the preced- 
ing by simply making a cyclical permutation of the subscripts attached to the G; 
firstly from x, y, z into y, z, x, secondly from x, y, z into z, x, y. On account of 
the divisor a 2 — 6 2 , small when the eccentricity of the planet's orbit is small, some 
of these formulas labor under a disadvantage. 

In computing the values of the definite integrals 

nab J p 3 ' nab J p s ' nab J p 3 
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where x, y, z are referred to the system of axes determined by those of the ellipse 
described by the attracting planet, we may suppose that the four quantities 
x , y , z , 7i are simply constants and unaffected by the transformation we have 
just made to pass from the system of coordinates x, y, % to that of x, Y, z. If, 
for the sake of discrimination, we designate the components along the mentioned 
axes as X', T 1 , Z', reserving X, T, Z for the components along the symmetric 
axes, we evidently have 

J' = al + a'F+ a"Z, T ~ (3X + 0' F+ fi'Z, Z' = yX + y'Y-\- y"Z. 

The components of the right members are determined the moment we have solved 
the cubic in G, for which the coordinates A, B, G furnish the requisite data. 
We have only to make 

x = A + ae, y Q = B, z = G, h = s/7P, 
where h is always to be taken positively. Thus 

X 1 - - a (A + ae) *L s i n 2 * — a'B J* cos 2 * + a" G ~ , 

Y' = -8(A + ae)M sin 2 * - B'B -Mr cos 2 * + p'cM , 

Z = -y(A + ae) ~~» sin 2 * — y>B -^-cos 2 * + y" O ^- . 

Substituting for a, a', a", B , B', B" their values in terms of y, y', y", if we put 

U x =-±±^ % M sin 2 *, 
* m* G 

n--~ ®jifL M cos 2 *, 

wr G 

we shall have yi — A f ^» 4- ^ _l ^ "1 

~ L G x + a 2 "*" G y + a 2 "*" Q, + a 2 J ' 

yi — B f ^ i Uy i ^ *] 

* -^L^ + J* ^ tf. + J 2 ^ <? 2 + &*J' 



Hill: Secular Perturbations of the Planets. 327 

By substituting in U x , U y , U^, the last values we have given for y, y', y", and 
bearing in mind that 

G z — G x = ™ 2 , G y — G x = m? sin 2 0, G z — G y = m 2 cos 8 0, 
the expressions for the first mentioned quantities become 

n _ A + ae . — • — txt-^ — r — r Ms'm z x 

u * ^F~ * G * & + a % G ' + h ) S iag ' 

U y = ~V-G y (G y + a»)( <?„ + &») M ™ % * X a , 
* m 5 » l ' ' A » ' 'sm^ cos 6 

U *= ^^^(^ + « a )(^ + 6 2 )^, 

where, if y, /, y" have been taken positively, the three radicals must receive the 
sign of G. 

Let it be required to find the component of this attracting force directed 
towards the centre of the ellipse described by the attracting planet. Putting 
r 2 = A* + B 2 + C*, we ought to multiply the components given above severally 

by the factors — ■ , — ■ , — , and take the sum. Which, if we do, and have 

regard to relations, G x , G y , G z satisfy as being the roots of the cubic in G; call- 
ing this component A, we have 

tA=U,+ U v + U.. 

But the component B directed towards the Sun will be more important. Put- 
ting r for the radius vector of the attracted planet, we have 

r\ = (A + aef + B 2 + G\ 
Then the components X, Y, Z ought to be multiplied severally by "*" ae , 

B G r ° 

— , , and thus is obtained 

rA=U x +U v +U z + aeX<. 

In fact, by multiplying X\ Y', Z' severally by the three systems of three mul- 
tipliers 

A + ae B G_ 

n n ' ~r^ 

B A + ae 







</rl—C 2 ' Vrl—G* ' 
(A + ae) G _ BG v/^— 77* 

''•„ Vr 2 - G 2 ' r VrT=~& ' ~~V ; 
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we shall arrive at a set of components, of which the first, already given, is 

directed toward the Sun, and the second is perpendicular thereto and lying in 

the plane of the attracting planet's orbit, while the third is perpendicular to both 

the preceding. If we call these new components X", 7", Z", and if the angle 

between the planes of the orbits of the two planets be denoted by J, we shall 

have 

R = X", S = 7" cos / — Z" sin /, W - 7" sin / + Z" cos J, 

where B is the component towards the Sun, S the component perpendicular 
thereto and lying in the plane of the attracted planet's orbit, and W the compo- 
nent perpendicular to that plane. But, if we make use of a rectangular set of 
elements instead of a polar, it is more commodious to refer the components to 
the plane and line of nodes of the attracted planet on the attracting planet's 
orbit. Let this ascending node be distant an arc = o from the perihelion of the 
latter. Then the desired components X'", 7'", Z'" will be 

X'" = X' cos a— F sino, 

T'" = X' cos /sin a + 7' cos I cos a — Z' sin /, 

Z'" = X' sin I sin a + T sin /cos a + Z' cos /. 

It is often interesting to know the position of the great circles forming the 
principal axes of the sphero-conic in reference to the great circle marked out in 
the heavens by the plane of the attracting planet's orbit. If we call Q, the lon- 
gitude of the ascending node of one of the planes of symmetry of the cone on 
the plane xy, and i the inclination (always between 0° and 180°), and r the 
angular distance of the centre of the sphero-conic from the node measured in 
the direction of increasing longitudes, the four following equations can be used 
for the determination of these quantities : 

tanisin ffi = -4 ^T^' ta * » cos Q = •£ a^T?> 






where the signs of the two radicals may be taken positive or negative, thus cor- 
responding to the four intei'sections of the two great circles with the one great 
circle. 
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To find the two semi-axes of the sphero-conic as measured by the arcs they 
subtend in the heavens, take the equation of the cone 

Lr x br y ir z 

the section of the cone by the plane z = V G z gives the ellipse whose equation is 

= 1. 



X 2 Y 2 



G r G, 



y 



The semi-axis of this in the direction of the axis, of x is V — G x and in the direc- 



tion of the axis of y, V — G v . The greatest latitude v] of the planet moving 
on the sphero-conic and the greatest longitude 2, will be given by the equations 



If I C 

tan>7 =y -^-, tanX = y — -^ . 



Or* T 0r s 

The general equation connecting the variables vj and "k will be obtained if, in the 
equation of the cone, we make 

x = p cos>7 sinX, Y = psin»7, z = p cos»7 cos A,, 
and thus is 

sin 2 % . tan 2 vj , cos 2 A, _ ft 
lr x ix y tr s 

These variables yj and % are expressed in terms of the perspective anomaly T as 
follows : 



cos y\ cos % = s V G z , cos sin A, = e V — G x cos T, sin yj = e V— G y sin T. 



/C 7t 
_2 • and when T= — , 

tan>? = y — -^- , and /l = 0; when T=n, >7 = 0,tana = — y-^-. 

On $e Solution of the Cubic in G and the Argument to be Employed in 

Tabulating the Elliptic Integrals. 

A few words may be added in reference to the solution of the cubic in G, 

* + *+£=!. 



G + a % ' G + b* ' G 
43 
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Using r 2 for A % -\- B 2 -\- C 2 , this equation expanded is 
£3 _ |y _ (« 2 + 6 2 )] G» + [a 2 ,! 2 + & 2 £ 2 - (a 2 + 6 s ) r 2 + « 2 & 2 ] — «W 2 = . 

It will be noticed that in obtaining the values of k, the modulus of the elliptic 
integrals involved, and of m,we do not need all the roots of this equation, but 
only their differences. Hence, it will be advantageous to put 

6?= J+ i [r 2 - (a 2 + 6 2 )] , a = a 2 + * [r 2 -(a 2 + 6 2 )] , 

b = 6 2 + \ [r 2 — (a 2 + & 2 )] , c = \ [r» - (a 2 + 6 2 )] 

The cubic will then take the form 



<7-j- a J -\-b J-\-c 

where /will serve us equally well as O ; but here we have 

a + b + c = A 2 + & + G\ 

and the developed equation in / will be 

J 3 + [a^L 2 + bB* + cC 2 — a 2 — b 2 — c 2 — ab — be — ca] / 

— [be A 2 + caJ5 2 + abC 2 — abc] = 0, 

By the elimination of c and G and the partial reintroduction of a 2 and b 2 , the 
shorter form is obtained, 

p _ [ a » + ab + b 2 — a? A 2 — b'B*] J+ a (¥& - ab) + b {a 2 A 2 — ab) = . 

Employing the well-known trigonometric process for the solution of the cubic 
having all its roots real, we derive ^ (between the limits =L 90°) from 

• o I — ^ a (Wff — ab) + b (a* A* - ab) 

The roots of the cubic are then 

J x - — 2 ry + ab + b 2 _ a s^ _ y^rj* gin (60° + ^) , 
v 3 

j y - -^[a' + ab + tf-afi'-^sinf 

j; = _L [a 2 + ab + b 2 — a % A 2 — WBJ sin (60° — •*) . 
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To render these equations simpler, we will put 
a + b = 2ra cos v, V 3 (a — b) = 2» sin v, a?A z + 6 2 -B 2 — 2 ab = nfy cos £, 

*= (a 2 J. 8 — 5 2 ^ 2 ) = n*n sin £. 



V 3 



With these modifications we have 



Sm3 ^ _ 2 (l-pooB*)l ! 



n(l — (UCOS^sin^O + 4-), 
» (1 — (i cos £)* sin 4* , 



V3 

2 
V3 



«£ = -^^(l —^ cos £ sin)* (60 — 4). 

The modulus of the elliptic integrals involved, k, will be given by the equation 
7{ s-„ in 2 e - J„ — J* = cos (60° - ^) 



«/« e/~ 



cosi^- 



and the quantity we have designated by m, by the formula 
m s = J z — /„ = 2n (1 — [i cos £) £ cos 4- 

The two elliptic integrals M sin 2 x, Jf cos 2 *: are functions of 4; consequently, 
they may be tabulated with the argument 4 or any function of 4 as, for instance, 
with sin 3^, that is, with the absolute discriminant of the cubic. We have 

2 = — i y— — ' , and, if we put x i for the second member of this, we have 

cos 4 



cos 



sin 



e ~] 



1 + */ cos 



= 4 



•05 



g + g 9 + g 85 + • • 



1 + x 1 + 2 (g 4 + g 16 + g 36 + ) 



The value of the nome g can be derived from the infinite series (three terms suffice 
except in very unusual cases), 

a = A 1 ~ x -4- — ( l — - V -4- — f * — g V _i_ 

2 * 1 + a; "'" 16 Vl -f a;/ ~ t ~ 512 Vl + xJ ' ' ' ' 
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As G x , G y , G z enter into many of the preceding formulas, it may be useful 
to note that 

G-* =J X + c, G y =J y + c, G t = J Z + c, 
G x + a? = J x + a, Gy + a z = J v +a, G, + a? = J t + a, 
G x + b* = J x + b, G y + b* = J v + b, G s + b* = J,+ b. 

Variation of the Elements of a Planet through Perturbation. 

The elements we select for use are defined thus : x, y, z, denoting the rectan- 
gular coordinates of a planet referred to the Sun as origin, adopt the elements 
<>x,Cy, c z ,f % ,f v ,f, such that 

ydz — zdy zdx — xdz xdy — ydx 

x ~ ~ dt ' ° v dt ' Cz dt ' 

/• o z dy — c y dz __ [ix * c x dz — o z dx (iy . c y dx — c x dy (iz 

Jx dt r ' Jv dt ~r~' Jz ~ "' ~dt T' 

where we note that p is the sum of the masses of the Sun and planet, and r is 
the radius vector of the latter. These six elements are not independent but 
satisfy the relation 

The additional element needed to complete the number six is the element every- 
where attached by addition to the time. With these constants the two equations 
of the path of the planet in space are 

c x x + c y y + c z z = 0, fir +f x x +f y y +f x z = cj + c\ -f- <% = &. 

To understand the correlation of the different terms of these equations it must 
be borne in mind that (i is a constant of three dimensions in reference to the 
linear unit. Consequently, the c are of two dimensions and the / of three dimen- 
sions in reference to this unit. 

The second equation belongs to a quadric surface not, in general, referred 
to its axes of symmetry. Removing the radical from it, it becomes 

(l/ -ft) x* + ({/ -/>) y> + {p* -/») 2? - 2/ x f y xy - 2f y f yz - 2/./. zx 

+ Wfe + 2¥f y y + 2Wfz — tt = 0. 
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The cubic, which must be solved in order that this may be referred to its axes 
of symmetry, is 

x s - [V- 07 +/v +ftl x* + W -ftW -ft + W -fl)W -ft 

But, if e denotes the eccentricity, we have/ 2 -f / 2 + f? = (i z e?; thus the preceding 
equation reduces to 

X 3 — ^(3 — ^)^ + ^(3 — 2e*) X — |« 6 (1— e 2 ) = 0. 

Put x — fay then 

x 3 — (3 — e 2 ) x 2 + (3 — 2e 2 ) x — (1 — e 2 ) = , 

or (* — l) 3 + e 2 (* — 1) 2 = 0. 

The three roots are x = 1, * = 1, * = 1 — e 2 . The second equation of the orbit 

of the planet is, therefore, a quadric of revolution about the major axis,, in 

the present investigation an ellipsoid, as we suppose e 2 <l. The expression 

c*+4 + c ^ s invariant, as also is /%+/$■{■/?, both being independent of the 

orientation of the axes of coordinates. The system of principal axes will be 

arrived at if we suppose f y = 0,/ z = 0, which imply also o x = 0, as is shown by 

the relation c x f x + c y f y + c z / z = 0, unless we have f x = Q, when, the quadric 

being a sphere, all systems of axes are principal. The equations of the orbit then 

take the form 

c y y + o z z = 0, fir + (*e% = & 2 , 

where & 2 = c 2 + c 2 . The radical removed, the second equation becomes 
l? (1 — e 2 ) x 2 + [i 2 y 2 + /mV + 2flfy«B — & 4 = 0. 

But, a being the semi-axis major, & 2 = ^a(l — e 2 ), therefore, the preceding 
equation takes the form 

I" + a -^4 +2 ^-<' -*> = •: 

or (te + ae) a , y* + g 2 __ ^ 

The adoption of the attracted planet instead of the centre of the Sun as the 
origin of coordinates renders it necessary, while employing the elements of the 
planet, to substitute — x , — y , — z for x, y, z. Thus, with this notation, the 
equations of the Sun's path in space are 

<Vi> + Wo + c > z o = > P r o — f* x o — fyVo - U* — #*• 
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The values of the differentials with respect to the time of the constants c 
and / must now be stated. Denoting the mass of the attracting planet by m', 

the perturbative function B for secular perturbations is — , and we have 
dc x dB dB m' , >. 

dt =yo d^- Zo w = T {Zoy ~ y ° z) ' 

dc„ dB dB m! , >. 

r 



dt "*dx ~"dz _ - 3 



dc z dB dB m! , . 

Tt =x »di -*sr = t<*0-"*)> 

dfz_dz l dcy__&yo de> _■ ™ ! (o .. „ „\ 
dt ~ dt dt dt dt "•" p 3 K * y yh 

df y dx dc z __ dz dc x , m! , ., 

~dT~~dT ~dt ~dT ~dT^"~f °* x) ' 
dt — dt dt dt dt ^y Kv h 

There is still one element to be added to the preceding to complete the number of 
six independent constants; this is the mean longitude at epoch. The well-known 
equation for its variation shows that a portion is immediately derivable from the 
motion of the perihelion, and another from that of the node ; what remains (call- 
ing the element T) is given by the equation 



_ an r dB . dB . dB~] m! r x , y , z -\ 



where n denotes the mean motion of the attracted planet. 

The system just given is the most general as respects the orientation of the 
axes of coordinates. Let us now specialize by taking, for the plane of xy, the 
plane of the orbit of the attracted planet. This makes o x = 0, c y = 0, c, = k 
and the equation between the c and the / shows that, in consequence, we have 
/, = 0. Then our equations become 

dc x _ m' df x _ dy dk , h m' 
~di~ ~J ~dt~ ~dt~dt ±k Y 
dc y m! dfy dx dk , m' 

dt ~ 9 s **' ~d7~ ~dt ~dt K y x ' 

dk _ w! i a df. __ dy dc x dx dc^ 

~dJ~ y ^ x ° yh It ~ ~dJ ~dJ ~dt ~dt' 

* For these formulas consult Laplace, Mecanique Celeste, Tom, I, Liv. II, Art. 64. 
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But our specialization, applied to the definitions of the /, gives 

These values enable us to -eliminate the differentials of x , y from the preceding 
equations. By differentiating the equation between the c and / we get 

» dc x . . dc y . , df z — 
f '~dt + /v ~dt +le ~dt-°- 

The substitution of preceding values in this renders it an identity, hence the 

df 
equation for -j± is superfluous. Moreover, after use in substitution, the equation 

dh 
for -=- is no longer needed, since the secular motion of the semi-axis a vanishes. 

Thus the equations to be employed are reduced to the five following: 

dc x , z 

w = -my oy> 



dc„ , z 

-&= mx *y 

%l~ ^L(f—n y±\(v x — x y \ — m'h x 

it- T\ /s ^ rj\ y " f x °yJ mk y 



If x denote the longitude of the perihelion of the attracted planet, we have 
/„ = fie cos x> fy — l^ sm % '■> a ^ s0 > ^ * ^ s the inclination and Q the longitude of 
the ascending node, o x = h sin i sin Q, c v = — h sin i cos Q ; whence 



d (sin i sin Q) __ ml z 

jt -k y *Y' 

d (sin i cos Q,) _ ml z 



dt k u p 3 

d(emn X ) — ^ f . y^y x_ y\_n^ k x 

— dt Tl esin r A y °p» Xo p 3 J y ?' 

d(eomv) m' / x \/ x y \ , tn' 7 y 

- ± -dt= nc-( e ™*-^A y °Y~*<-f) + T h Y-- 



336 Hill : Secular Perturbations of the Planets. 

The point of departure for longitudes is still undetermined ; but, if this is 
taken to be the ascending node of the attracted planet's orbit on that of the attract- 
ing planet, we can employ the X'", T'", Z"' we have previously denned. Thus, using 
brackets with subscript to denote integration along the orbit of the attracting 
planet we have the equations 



[ 



•d(s'mi sin£)-| _ ml ym 

- ~ ~dt Jo ¥ y ° z ' 

^ , — _ JUL. sr„ Z"> 

dt Jo h ° ' 



/// 



in 



d(sini cos &) n _ __ m' 
dt Jo h 

[^1 = 4(.-.*--fe-)(* i --^)-T-*- r 

["-^1= 4-C— *- J £X**"-+ rm ) + T kYI 

The integration round the orbit of the disturbed planet has still to be exe- 
cuted in order to arrive at the secular motion of the elements. For this we are 
confined to the use of mechanical quadratures. Here we may use either of the 
three anomalies, or any variable which will show the position of the planet on 
its orbit, as the independent variable. 
West Nyack, March 13, 1901. 



